Five-loop renormalization-group expansions 
for the three-dimensional n-vector cubic model 
and critical exponents for impure Ising systems 

O" 

o 

(N) ■ D. V. Pakhnin 1 and A. I. Sokolov 1 ' 2 

1 Department of Physical Electronics, Saint Petersburg Electrotechnical University, 
Professor Popov Street 5, Saint Petersburg 197376, Russia, 
If} \ 2 Department of Physics, Saint Petersburg Electrotechnical University, 

Professor Popov Street 5, Saint Petersburg 197376, Russia 

O 

S ■ Abstract 

> . 

^ ■ The renormalization-group (RG) functions for the three-dimensional n- 

vector cubic model are calculated in the five-loop approximation. High- 
precision numerical estimates for the asymptotic critical exponents of the 
three-dimensional impure Ising systems are extracted from the five-loop 
£h ■ RG series by means of the Padc-Borcl-Lcroy rcsummation under n = 0. 

These exponents are found to be: 7 = 1.325 ± 0.003, r] = 0.025 ± 0.01, 
v = 0.671 ± 0.005, a = -0.0125 ± 0.008, and (3 = 0.344 ± 0.006. For the 
\ correction-to-scaling exponent, the less accurate estimate u = 0.32 ± 0.06 is 
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I. INTRODUCTION 



The critical thermodynamics of cubic crystals and weakly disordered systems has re- 
mained an area of extensive theoretical work during past decades. Considerable progress in 
studying the random critical behavior was achieved 25 years ago when Harris and Lubensky 
|T|fJ] and Khmelnitskii [H] attacked the problem by the field-theoretical renormalizat ion- 
group (RG) approach based on the Euclidean scalar (p A theory in (4 — e) dimensions. As 
a results, the regular method for calculating critical exponents and other universal quanti- 
ties of the impure Ising model - the famous ^-expansion - was invented. The numerical 
power of this technique, however, stayed for a long time unclear since only lower-order con- 
tributions to critical exponents and the equation of state have been found 0-0]. Recently, 
starting from the five-loop RG series obtained for the (4 — e)-dimensional cubic model by 
Kleinert and Schulte-Frohlinde ||, the calculation of the A/e-expansions for critical expo- 
nents was performed up to the y/e 4 and \fe terms |]]. As was found, these series possess a 
rather irregular structure making them unsuitable for subsequent resummation and, hence, 
practically useless for getting numerical estimates |J. 

On the other hand, there exists an alternative field-theoretical approach that proved 
to be very efficient when used for evaluation of universal critical quantities. We mean 
the perturbative renormalization group in three dimensions (3D) yielding most accurate 
numerical estimates for critical exponents, critical amplitude ratios, and universal higher- 
order couplings of the 0(n)-symmetric systems P-pOf. The impure Ising model at criticality 
is known to be described by the n- vector field theory with the quartic self-interaction having 
a hupercubic symmetry, provided n — > (the replica limit) and the coupling constants 
have proper signs. In the 1980s, the RG expansions for three-dimensional (3D) cubic and 
impure Ising models have been calculated in the two-loop |2~1 |, three-loop p2| , |2~3|j , and four- 
loop p4| , p5l approximations paving the way for estimating the critical exponents and other 
universal quantities fl2~l-|34|j. The numerical results thus obtained was found to agree, in 
general, with the most accurate experimental and simulation data. 

In the course of this study, it was revealed, however, that even the highest-order available, 
four-loop 3D RG expansions, when resummed by means of the generalized Pade-Borel-Leroy 
method do not allow us, in fact, to optimize the resummation procedure, i. e., to choose the 
best Pade approximant and the optimal value of the tune parameter since there is the only 
approximant - [3/1] - that does not suffer from positive axis poles. Moreover, an account 
for four-loop terms in the 3D RG expansions shifts the fixed point coordinates and the value 
of the correction-to-scaling exponent u appreciably with respect to their three-loop analogs, 
indicating that at this step the RG-based iterations do not still achieve their asymptote. 
This prevents the four-loop RG approximation from to be thought of as sufficient, i. e., 
providing, within the perturbation theory, the accurate theoretical predictions. 

In such a situation, a calculation of the higher-order contributions to the RG functions 
looks very desirable. In this paper, the five-loop RG expansions for the three-dimensional 
cubic model are obtained and resulting numerical estimates for the critical exponents of the 
weakly disordered Ising systems are found. 
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II. RG EXPANSIONS FOR FUNCTIONS AND CRITICAL EXPONENTS 



The Landau- Wilson Hamiltonian of the three-dimensional n-vector cubic model reads: 
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where <^ is an n-component real order parameter, m 2 being the reduced deviation from 
the mean-field transition temperature. In the replica limit, this Hamiltonian describes the 
critical behavior of the impure Ising model provided u Q < and v Q > 0. 

We calculate the (3 functions for the Hamiltonian Eq. ([J) within a massive the- 
ory. The renormalized Green function Gn(p,m), the if 2 insertion and four-point vertices 
Ur(pz, m, u,v,), Vn(pi, m, u, v, ) are normalized at zero external momenta in a conventional 
way: 



G^O.m) = m 2 , 
r ( R 2 \p,q,m,u,v) 



dG R \p,m) 
dp 2 



p 2 =0 



1. 



p=g=0 



Un(0,m,u,v) — mu, Vr(0, m, u, v) — mv . 



(2) 



The value of the one-loop vertex graph including the factor (n + 8) is absorbed in u and v 
in order to make the coefficient for the u 2 term in (3 U equal to unity. 



The four-loop RG expansions for the functions of interest have been found earlier 24 



To extend these series to the five-loop order, we calculate corresponding tensor (field) factors 
generated by the 0(n)-symmetric and cubic interactions. Taking, then, numerical values of 
the 3D integrals from Ref. EH], we arrive to the following five- loop expansions: 
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These expansions will be used to evaluate the critical exponents of the impure Ising model. 



III. RESUMMATION AND NUMERICAL ESTIMATES 



Numerical values of critical exponents are known to be determined by the coordinates 
of a relevant fixed point. In our case, i. e., for n = 0, a point of interest is the random fixed 
point. To find its location and, more generally, to extract the physical information from 
divergent RG series, a proper resummation procedure should be applied. Here, we use the 
Pade-Borel-Leroy resummation technique, which demonstrates high numerical effectiveness 
both for the 0(n)-symmetric models |P| JTT| , [T5[] and for anisotropic systems preserving their 
internal symmetries (see, e. g. Ref. f3"6[| ) for detail). Since the expansions of quantities 
depending on two variables u and v are dealt with, the Borel-Leroy transformation is taken 
in a generalized form: 

To perform an analytical continuation, the resolvent series 

is constructed, which is a series in powers of A with coefficients being uniform polynomials 
in x, y and then Pade approximants [L/M] in A at A = 1 are used. 

For the resummation of the five-loop RG expansions, we employ three different Pade 
approximants: [4/1], [3/2], and [2/3]. The first of them, being pole-free, is known to give 
good numerical results for 3D 0(n)-symmetric models while the others are near-diagonal and 
should reveal, a priory, the best approximating properties. The coordinates of the random 
fixed point resulting from the series Eqs. (|^, ||) resummed using these approximants under 
6 = and 6=1 are presented in Table I, which also contains analogous estimates given by 
the four-loop RG expansions. The four-loop series were processed on the base of the Pade 
approximant [3/1], since use of the diagonal approximant [2/2] leads to the integrand in 
Eq. (^) that has a dangerous pole in the vicinity of the random fixed point both for (3 U and 
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(3 V |37| . The fixed-point location given by the approximant [2/3] is presented for b = only, 



because for 6=1 this approximation predicts no random fixed point. 

As is seen from Table I, Pade approximants [4/1] and [3/2] yield numerical values of the 
random fixed-point coordinates, which are remarkably close to each other. Moreover, for b = 
they are also close to those given by the approximant [3/1]: the largest difference between 
the five-loop and four-loop estimates does not exceed 0.026. With increasing b, corresponding 
numbers diverge, indicating that b = is an optimal value of the tune parameter. On 
the contrary, Pade approximant [2/3] gives the random fixed-point location, which deviates 
appreciably from those predicted by approximants [4/1], [3/2], and [3/1]. This approximant, 
however, is found to lead to poor numerical results even for simpler systems, e. g., for the 
Ising model. Indeed, when used to evaluate the coordinate of the Ising fixed point, it results 



in v c = 1.475 (under b = 0) while the best estimate today is known to be v c = 1.411 jjTJ 
This forces us to reject the data obtained on the base of the approximant [2/3]. 

So, to determine the random fixed-point coordinates, we have to average the numerical 
data given by three working Pade approximants at b = 0. This procedure yields the values 

u c = -0.711, v c = 2.008, (9) 

which are claimed to be the final results of our search of the random fixed-point location. 
Let us estimate their accuracy. It seems unlikely that the deviations of these numbers from 
the exact ones would exceed the differences between them and the four-loop estimates since, 
among all proper estimates, the four-loop ones most strongly differ from the averaged values 
Hence, the error bounds for u c and v c are believed to be not greater than ±0.012 and 
±0.016, respectively. Another way to estimate an apparent accuracy is to trace how the 
averaged values of the random fixed point coordinates vary with the variation of b. We 
calculate u c and v c using the pole-free approximants [4/1] and [3/1] for < b < 15; b = 15 is 
chosen as a largest reasonable value of the tune parameter, since for greater b the saturation 
of dependences of various quantities on b becomes visible. Running through this interval, 
the averaged coordinates change their values by about 0.02 indicating that an accuracy of 
the estimates (|9]) is of the order of 0.01-0.02, in accord with that found above. 

With the numbers ([5]) in hand, we can evaluate the critical exponents for the 3D impure 
Ising model. The exponent 7 is estimated by the Pade-Borel-Leroy resummation of the 
RG series (|5|) for 7 _1 and of the inverse series, i. e., the RG expansion for 7. The Fisher 
exponent is also evaluated in two different ways: via the estimation of the critical exponent 
772 = (2 — 77) ('y -1 — 1) having the RG expansion, which exhibits a good summability, and by 
direct substitution of the fixed point coordinates into the series (§) with rapidly diminishing 
coefficients. The estimates for rj originating from 772 were obtained under the central value 
of 7: 7 = 1.325. The numerical results thus found are collected in Table II. 

As is seen from this Table, two methods of evaluating of the susceptibility exponent 7 
lead to remarkably close numerical results, which very weakly depend on the tune parameter. 
Indeed, with increasing b from to 15, the estimates for 7 obtained by the resummation 
of the RG series for 7 and 7 _1 on the base of the pole-free approximant [4/1] vary by 
less than 0.0036, while the difference between them never exceeds 0.0013. Under the same 
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variation of b, the value of 7 averaged over these two most reliable approximations remains 
within the segment [1.3240, 1.3266]. On the other hand, the accuracy of determination 
of the critical exponents depends not only on a quality of the resummation procedure but 
also on the accuracy achieved in the course of locating of the relevant fixed point. That is 
why we investigated to what extent the numerical estimates for 7 vary when coordinates of 
the random fixed point run through their error bars. It was found that the susceptibility 
exponent calculated at the optimal value of tune parameter 6 = 2 (see Table II) does not 
leave the segment [1.3228, 1.3263]. This enables us to conclude that the error bounds for 
the theoretical value of 7 obtained in this work would be about ±0.003 or smaller. 

Less stable numerical results, with respect to a variation of b, are found for the Fisher 
exponent 77. As one can see from Table II, the values of 77 given by the resummation of the 
RG series for 772 with use of the pole-free Pade approximant [4/1] spread from 0.0148 to 
0.0312. The average over this interval is equal to 0.023. The direct summation of the RG 
expansion for 77 at the random fixed point gives 0.027. It is natural to suppose that 0.025 
should play a role of the most likely value of exponent rj. Since the estimates for 77 found via 
the evaluation of 772 are sensitive to the accepted value of 7, the apparent accuracy achieved 
in this case is not believed to be better than ±0.01. 

Having estimated 7 and 77, we can evaluate other critical exponents using well-known 
scaling relations. So, the final results of our five-loop RG analysis are as follows: 

7 = 1.325 ±0.003, 77 = 0.025 ±0.01, v = 0.671 ± 0.005, 

a = -0.0125 ±0.008, (3 = 0.344 ± 0.006. (10) 

These numbers are thought to be the most accurate theoretical estimates for the critical 
exponents of the 3D impure Ising model known today. It is interesting to compare them 
with those obtained earlier within the lower-order RG approximations. For the exponent 7, 



previous RG calculations in three dimensions gave the values 1.337 pl|j27]] (two-loop), 1.328 



(three-loop), 1.326 [|4j] (four-loop), and 1.321 |25|] (four- loop). Being found by means of 



the different resummation procedures, they are, nevertheless, centered around our estimate 
which is thus argued to be very close to the exact value of 7 or, more precisely, to the true 
asymptote of the RG iterations. 

At the end of this work, we employ our technique to evaluate the correction-to-scaling 
exponent uo. The exponent u is known to be equal to the stability matrix eigenvalue that 
has a minimal modulus. The derivatives df3 u /du, dfi u /dv, d(3 v /du, and df3 v /dv entering this 
matrix are evaluated numerically at the random fixed point on the base of the resummed 
RG expansions for j3 u and (3 V , and then the matrix eigenvalues are found. Such a procedure 
leads to the estimates for uj presented in Table I (lower lines). They are seen to be consid- 
erably scattered and sensitive to the tune parameter. The average over three working Pade 
approximants, however, being equal to 0.315 at b = and to 0.316 at b = 1 turns out to be 
stable under the variation of b unless b becomes large. It is natural therefore to accept that 

to = 0.32 ±0.06. (11) 

This number is smaller by 0.05 — 0.07 than its counterparts given by recent Monte Carlo 
simulations [B8J and the alternative RG analysis pEl, but their central values lie within the 
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declared error bounds (|XI|) . Hence, an agreement between the results discussed exists. On 
the other hand, because of its low accuracy, the estimate ([TlD would not be thought of 
as satisfactory. It certainly needs to be improved, along with the estimates for the small 
exponents rj and a also exhibiting appreciable uncertainties. The only way to do this is 
the studying of the critical behavior of the 3D impure Ising model in the next perturbative 
order that includes calculations of the six-loop terms in the relevant RG expansions. Such 
calculations are now in progress. 



IV. CONCLUSION 

To summarize, we have calculated the five-loop RG expansions for the j3 functions and 
critical exponents of the 3D n-vector cubic model. The resummation of the RG series 
by the Pade-Borel-Leroy technique in the replica limit (n = 0) has enabled us to obtain 
high-precision numerical estimates for the random fixed-point coordinates and the critical 
exponents 7, rj, a, u, and (3 of the 3D impure Ising model. For the correstion-to-scaling 
exponent u>, the resummed five-loop RG expansions turned out to give much less accurate 
numerical results. The values of the critical exponents obtained earlier from the lower-order 
RG expansions in three dimensions are shown to be centered by our five-loop estimates. 
This indicates that the five-loop RG approximation provides numerical data very close to 
the asymptotic ones, i. e., to those representing the point of convergence of the RG-based it- 
erations. At the same time, an accuracy achieved when evaluating the exponent u would not 
be referred to as sufficient making the next-order, six-loop RG calculations, very desirable. 

Note added. After this paper had been submitted for publication, Ref. |[{| appeared 



where the six-loop RG expansions for the 3D cubic model are calculated. The five-loop terms 
obtained in this remarkable work are found to agree with ours. Among other quantities, 



in Ref. |39| the marginal value of n, n c , separating different regimes of critical behavior of 
the cibic model, is estimated. To clear up how sensitive to the method of resummation 
the five-loop RG results are, we evaluate n c using Pade-Borel-Leroy technique. Exploiting 
the near-diagonal Pade approximant [3/2] under b varying from to 20, the values of n c 
are obtained which lie between 2.89 and 2.92. They agree quite well with the estimate 
n c = 2.91(3) extracted from the five-loop RG series by means of the conformal-mapping- 
based resummation machinery [|3£ . 
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TABLES 



TABLE I. Numerical estimates for the random fixed-point location and correction-to-scaling 
exponent uj obtained from the five-loop RG expansions |3|) resummed by the Pade-Borel-Leroy 
technique using approximants [4/1], [3/2], and [2/3]. The last column contains results given by the 
four-loop RG series processed on the base of the approximant [3/1]. The superscript "c" denotes 
that the exponent uj is complex and its real part is presented. The superscript "p" stands to mark 
that the Pade approximant has a " nondangerous" positive axis pole, i. e. a pole well remoted from 
the origin (t > 40) that affects neither the procedure of numerical evaluation nor the value itself 
of the Borel integral. 





b 


[4/1] 


[3/2] 


[2/3] 


[3/1] 







-0.7200 


-0.7148 


-0.6871 


-0.6991 




1 


-0.7445 


-0.7385? 




-0.6839 


V c 





2.0182 


2.0125 


2.0571 


1.9922 




1 


2.0296 


2.0236? 




1.9877 


it! 





0.266 


0.303 


0.462 c 


0.376 




1 


0.263 


0.325? 




0.361 



TABLE II. Numerical estimates for the critical exponents 7 and 77 obtained from the five-loop 
RG expansions (||, ||) resummed by the Pade-Borel-Leroy technique using approximants [4/1] and 
[3/2]. DS stands for "direct summation", the symbol (7~ 1 ) -1 means that the RG series for 7 _1 
was resummed. The superscript "p" denotes, as in Table 1, that the Pade approximant has a 
"nondangerous" pole, while empty cells are due to the dangerous ones spoiling corresponding 
approximations. The estimates for 77 standing in the fifth and sixth lines were obtained under 
7 = 1.325 by the resummation of the RG series for 772. 



b 







1 


2 


3 


4 


5 


10 


15 




[4/1] 


1.3236 


1.3244 


1.3250 


1.3254 


1.3257 


1.3260 


1.3268 


1.3272 




[3/2] 








1.3253? 


1.3257? 


1.3260 


1.3265 


1.3267 


7 


[4/1] 


1.3245 


1.3248 


1.3250 


1.3252 


1.3253 


1.3254 


1.3257 


1.3259 




[3/2] 


1.3246? 


1.3251? 


1.3254? 


1.3257? 


1.3259? 


1.3261? 


1.3267? 


1.3270? 


77 (via 7/2) 


[4/1] 


0.03122 


0.02765 


0.02506 


0.02311 


0.02158 


0.02036 


0.01665 


0.01478 




[3/2] 








0.02870? 


0.02419? 


0.02168? 


0.01666? 


0.01487 


r, (DS) 










0.0272 
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